REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  the  burden,  to  the  Department  of  Defense,  Executive  Service  Directorate  (0704-0188).  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no 
person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ORGANIZATION. 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORTTYPE 

30-09-2011  Final 


4.  TITLE  AND  SUBTITLE 

MULTI-SCALE  PREDICTION  AND  SIMULATION  OF  LOCALIZATION 
BANDING  IN  GRANULAR  MEDIA 


3.  DATES  COVERED  (From  -  To) 

15/04/08-30/06/11 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

FA9550-08-1-0192 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Jose  E.  Andrade 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Northwestern  University 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Office  of  Scientific  Research 
875  North  Randolph  Street 
Suite  325,  Rm  31 12 
Arlington,  VA  22203 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

AFOSR 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

AFRL-OSR- VA-TR-201 2-0254 


14.  ABSTRACT 

Shear  bands  are  a  critical  factor  in  the  response  of  granular  materials,  such  as  soils,  rocks,  and  concrete.  The  objective  of  this  research  is  to  develop 
methods  for  accurately  predicting  the  onset  and  propagation  of  deformation  bands  in  granular  media  under  extreme  environments  by  developing  an 
efficient  multiscale  simulation  framework.  The  goal  is  to  bridge  continuum  scale  methods  with  the  existing  discrete  methods  to  obtain  an  accurate 
multiscale  model  that  can  be  used  to  predict  and  model  shear  band  inception  and  propagation  in  granular  media.  Fundamental  understanding  of 
granular  mechanics  under  shear  banding  will  pave  the  way  to  tackling  more  challenging  phenomena  such  as  projectile  penetration  and  blasting. 


15.  SUBJECT  TERMS 

Penetration,  granular  matter,  multiscale  modeling,  numerical  methods 


16.  SECURITY  CLASSIFICATION  OF:  1 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

u 

u 

u 

17.  LIMITATION  OF  18.  NUMBER  |19a.  NAME  OF  RESPONSIBLE  PERSON 


ABSTRACT 


OF 

PAGES 


|19b.  TELEPHONE  NUMBER  (Includearea  code) 


Standard  Form  298  (Rev.  8/98) 

Prescribed  byANSI  Std.  Z39.18 

Adobe  Professional  7.0 


MULTISCALE  PREDICTION  AND  SIMULATION  OF  LOCALIZATION  BANDING 

IN  GRANULAR  MEDIA 


FA9550-08-1-0192 


Jose  E.  Andrade 

Division  of  Engineering  and  Applied  Science 
California  Institute  of  Technology,  CA 

Ted  Belytschko 

Department  of  Mechanical  Engineering 
Northwestern  University,  Evanston,  IL 


Abstract 

Shear  bands  are  a  critical  factor  in  the  response  of  granular  materials,  such  as  soils,  rocks, 
and  concrete.  The  objective  of  this  research  is  to  develop  methods  for  accurately 
predicting  the  onset  and  propagation  of  defonnation  bands  in  granular  media  under 
extreme  environments  by  developing  an  efficient  multiscale  simulation  framework.  The 
goal  is  to  bridge  continuum  scale  methods  with  the  existing  discrete  methods  to  obtain  an 
accurate  multiscale  model  that  can  be  used  to  predict  and  model  shear  band  inception  and 
propagation  in  granular  media.  Fundamental  understanding  of  granular  mechanics  under 
shear  banding  will  pave  the  way  to  tackling  more  challenging  phenomena  such  as 
projectile  penetration  and  blasting. 

Status/Progress 

A  key  component  for  the  success  of  the  multiscale  scheme  is  the  ability  to  extract  the 
continuum  behavior  directly  from  the  granular  scale.  To  achieve  this  goal,  a  number  of 
objectives  have  been  successfully  achieved:  1.  The  design  and  implementation  of  a  semi¬ 
concurrent  multiscale  scheme  to  extract  material  behavior  from  the  granular  scale  and 
pass  it  onto  the  continuum  scale  [1],  2.  The  creation  of  a  novel  semi-implicit  return 
mapping  algorithm  to  enable  the  semi-concurrent  multiscale  computations  [2],  3. 
Development  of  constitutive  models  capable  of  capturing  rate  effects  [3],  4.  Design  and 
implementation  of  parallel  algorithms  to  perform  multi-unit  cell  calculations  [4,  6],  In 
addition  to  these  achievements,  final  efforts  focused  on  crucial  aspects  dealing  with  the 
development  of  a  new  hierarchical  multiscale  method  amenable  to  direct  linkage  with 
experimental  and  numerical  granular  kinematics  [6],  and  the  development  of  multiscale 
failure  modeling  methods  linking  assumed  enhanced  strain  (AES)  finite  elements  [5]  with 
discrete  element  calculations.  In  this  final  report  for  this  project,  we  focus  on  the  overall 
achievements. 

Perhaps  the  two  most  important  achievements  in  this  project  were:  1)  development  of  the 
first  predictive  (validated)  multiscale  method  for  failure  modeling  in  granular  materials 
and  2)  multiscale  enhancement  of  existing  failure  modeling  techniques  such  as  Assumed 
Enhanced  Strains  (AES)  to  accurately  model  multiscale  failure  in  granular  solids.  These 


achievements  enabled  us  to  accomplish  the  goals  proposed  for  this  project,  which  are 
shown  pictorially  in  Figure  1. 


Figure  1.  Schematic  of  proposed  research  and  its  main  goal:  to  enhance  elements  with 
embedded  discontinuities  within  a  multiscale  framework  connecting  continuum  and 
granular  scales. 


To  showcase  the  aforementioned  achievements,  we  will  focus  on  the  last  publication 
conceived  during  this  project  [5].  In  this  paper,  we  developed  a  new  Assumed  Enhanced 
Strain  (AES)  technique,  using  an  original  formulation  proposed  in  [7].  The  original 
formulation  enriched  elements  but  could  not  describe  plasticity  inside  shear  bands 
accurately.  Our  new  method  accesses  the  granular  scale  to  enhance  the  constitutive 
description  at  the  AES  continuum  level.  Figure  2  shows  an  example  of  simulations 
performed  using  discrete  mechanics  and  compared  with  simulations  using  the  AES 
multiscale  technique. 


Figure  2.  Sample  (left)  composed  of  regularly  packed  granular  materials  under  plane 
strain  compression.  A  shear  band  develops  when  vertical  strain  is  about  0.5%,  with  this 
shear  band  dominating  the  process  for  the  duration  of  the  discrete  element  simulation. 


We  used  our  proposed  AES  technique  to  simulate  the  aforementioned  sample  of  granular 
material  under  the  same  boundary  conditions.  Figure  3  shows  the  results  using  a 
triangular  irregular  mesh  enhanced  with  AES  elements  that  have  multiscale  capabilities 
to  describe  the  constitutive  response  based  on  the  grain  scale  mechanics.  This  figure  also 
shows  a  comparison  between  the  global  stress-strain  responses  between  the  multsicale 
AES  and  the  direct  numerical  simulation  using  discrete  mechanics  (DEM).  Clearly,  the 
method  is  able  to  capture  the  most  salient  features  related  to  overall  strength  of  the 
sample  and  the  topology  of  the  shear  band. 
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Figure  3.  Results  comparing  multiscale  AES  with  direct  simulations  using  discrete 
mechanics  (DEM). 

Finally,  we  have  shown  that  the  method  can  be  validated  using  actual  experiments.  This 
is  the  first  method  of  its  kind  to  be  ever  compared  with  actual  physical  experiments.  In 
reference  [5],  we  compare  the  multiscale  AES  method  with  physical  experiments  on 
sands  where  a  dominant  shear  band  was  observed.  Micromechanical  information 
measured  during  the  test  was  used  to  construct  the  hierarchical  multiscale  AES 
calculation.  Figure  4  shows  the  results  of  the  multiscale  AES  prediction  compared  with 
experimental  data.  It  can  be  seen  that  the  method  clearly  captures  the  macroscopic 
response  of  the  sample  accurately. 


Figure  4.  Comparisons  of  multiscale  AES  with  physical  experiments  where  a  dominant 
shear  band  was  observed. 
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